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Abstract 


We derive extremal black hole solutions for a variety of four dimensional models 
which, after Kaluza-Klein reduction, admit a description in terms of 3D gravity coupled 
to a sigma model with symmetric target space. The solutions are in correspondence 
with certain nilpotent generators of the isometry group. In particular, we provide the 
exact solution for a non-BPS black hole with generic charges and asymptotic moduli 
in NV = 2 supergravity coupled to one vector multiplet. Multi-centered solutions can 
also be generated with this technique. It is shown that the non-supersymmetric solu- 
tions lack the intricate moduli space of bound configurations that are typical of the 
supersymmetric case. 
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1 Introduction 


Soon after the attractor mechanism was first discovered in supersymmetric (BPS) black holes 
[1], it was reformulated in terms of motion on an effective potential for the moduli [2]. Fer- 
rara et al demonstrated that the critical points of this potential correspond to the attractor 
values of the moduli. More recently, several groups used the effective potential to show that 
non-supersymmetric (non-BPS) extremal black holes can also exhibit the attractor mecha- 
nism, thereby creating a new and exciting field of research [3] 4]. Many connections between 
non-BPS attractors and other active areas of string theory soon revealed themselves. An- 
drianopoli et al found that both BPS and non-BPS black holes embedded in a supergravity 
with a symmetric moduli space can be studied using the same formalism, and they uncov- 
ered many intricate relations between the two [5] [6]. Dabholkar, Sen and Trivedi proposed 
a microstate counting for non-BPS black holes (albeit subject to certain constraints [7]). 
Saraikin and Vafa suggested that a new extension of topological string theory generalizes 
the Ooguri-Strominger-Vafa (OSV) formula such that it is also valid for non-supersymmetric 
black holes [8]. Studying non-BPS attractors could also give insight into non-supersymmetric 
flux vacua. Given all these possible applications, it is important to characterize non-BPS 
black holes as fully as possible. 


There has been a great deal of progress in understanding the near-horizon region of these 
non-BPS attractors. The second derivative of the effective potential at the critical point 
determines whether the black hole is an attractor, and the location of the critical point 
yields the values of the moduli at the horizon; in this way, one can compute the stability 
and attractor moduli for all models with cubic prepotential [9] [10]. However, the effective 
potential has only been formulated for the leading-order terms in the supergravity lagrangian. 
If one wants to include higher-derivative corrections, one can instead use Sen’s entropy 
formalism, which incorporates Wald’s formula, to characterize the near-horizon geometry in 
greater generality [II]. Sen’s method has led to many new results [12 {13} [14] (15) [16] [17]. 
The tradeoff is that this method cannot be used to determine any properties of the solution 
away from the horizon. 


The BPS attractor flow is constructed from the attractor value z%ps = z%ps(p", qr) by 
simply replacing the D-brane charges with the corresponding harmonic functions: 


Zpps(Z) = Zpps(P* = H' (2), qr > H;(2)) (1.1) 
where the harmonic functions are 
H! (T) h? 1 p! ) 
> = + Toa. 
H(z) hy |Z] \ ar 


Moreover, this procedure can be applied to construct the multi-centered BPS attractor flow 
that describes the supersymmetric black hole bound state , where the harmonic functions 


are generalized to have multiple centers: 
H' (z h! 1 D; 
( (2) \ _ +a = (1 
A; (#) hy ~ea \ (ar) 
It is conjectured in that the non-BPS flow can be generated in the same fashion, namely, 
by replacing the charges in the attractor value with the corresponding harmonic functions. 


However, as will be proven in this paper, this procedure does not work for systems with 
generic charge and asymptotic modulit] 


In principle, one could construct the full non-BPS flow (the black hole metric, together 
with the attractor flow of the moduli) by solving the equation of motion derived from the 
lagrangian. However, this is a second-order differential equation and only reduces to a first- 
order equation upon demanding the preservation of supersymmetry. Ceresole et al have 
written down an equivalent first-order equation in terms of a “fake superpotential”, but 
so far, the fake superpotential can only be explicitly constructed for special charges and 
asymptotic moduli 20]. The most generic non-BPS equation of motion is complicated 
enough that it has not yet been solved. Similarly, multi-centered non-BPS black holes have 
not been studied. 


Our goal is to construct the full flow for non-BPS stationary black holes in four di- 
mensions. Instead of directly solving the equation of motion, we reduce the action on the 
timelike isometry and dualize all 4D vectors to scalars. The new moduli space M3p contains 
isometries corresponding to all the charges of the black hole, and the black hole solutions 
are simply geodesics on M3p. This method was introduced in and has been used to 
construct static and rotating black holes in heterotic string theory and to study the 
classical BPS single-centered flow and its radial quantization [26]. 


In this paper, we work in two specific theories of gravity, but we expect that this method 
can be used for any model whose Msp is symmetric. The basic technique is reviewed in 
more detail in Section 2. In Section 3, we show how this method works in a simple case: the 
toroidal compactification of D-dimensional pure gravity. Section 4 serves as an introduction 
to single-centered attractor flow in M = 2 supergravity coupled to one vector multiplet, and 
Sections 5 and 6 are dedicated to constructing the full flows for both BPS and non-BPS 
single-centered black holes with generic charges. We find that they are generated by the 
action of different classes of nilpotent elements in the coset algebra. Both types of flows are 
shown to reach the correct attractor values at the horizon. In Section 7, the procedure is 
generalized straightforwardly to construct both BPS and non-BPS multi-centered solutions. 
We use a metric ansatz with a flat spatial slice and we are able to recover the BPS bound 
states described by Bates and Denef. Using the same ansatz, we are able to build non-BPS 
multi-centered solutions. Unfortunately, solutions generated this way turn out to always have 


'This has also been shown in [20]. 


charges at each center which are mutually local. The last section reviews our conclusions 
and suggests possibilities for future work. 


2 Framework 


Here we outline the method we will use to construct black hole solutions. This method was 
first described in [22]. We first reduce a general gravity action from four dimensions down to 
three, and derive the equation of motion. We then specialize to certain theories which have 
a 3d description in terms of a symmetric coset space. In such situations, we can easily find 
solutions to the equation of motion. The solutions are geodesics (or generalizations thereof) 
on the 3d moduli space and they are generated by elements of the coset algebra. 


2.1 3D Moduli Space 


We will study stationary solutions in a theory with gravity coupled to scalar and vector 
matter. Let the scalars be zë and the vectors be A’. Then the most general ansatz for a 
stationary solution in four dimensions is: 


ds? —e” (dt +w)? + e g dr"dr® (2.1) 
F! = dA'=d(Aj(dt+w) +A‘) 


where a,b = 1, 2,3 label the spatial directions and bold font denotes three-dimensional fields 
and operators. Since none of the fields are time-dependent, we can compactify on the time 
isometry and reduce to three-dimensional space M3p. This procedure is called the c*-map. 
In three dimensions, a vector is Hodge-dual to a scalar. The equations of motion for w and 
the gauge fields allow us to define the dual scalars ¢,, and @ar. 


We then obtain the 3d lagrangian in terms of only scalars 


1 1 
L= zv8(-7zR + Oa OH" Jmn) (2.3) 
where ġ” are the moduli fields 
g” = {U, 21,2, bu Ads Par} (2.4) 


and ga is the space time metric and gmn is the metric of a manifold Msp. The system is 
3D gravity minimally coupled to a nonlinear sigma model with moduli space Msp. Next, 
we will find the equation of motion in this theory. 


2.2 Attractor Flow Equation 


The equation of motion of 3D gravity is Einstein’s equation: 


1 1 
Rap — z8ak = KTab = K(O.0" O60" Jmn a Babe" OO" Jmn) (2.5) 


and the equation of motion of the moduli is: 
VaVe” + T hapla” = 0 (2.6) 


For simplicity, we consider only the case where the 3D spatial slice is flat (it is guaranteed 
to be flat only for extremal single-centered black holes). Then the dynamics of the moduli 
are decoupled from that of the 3D gravity: 


Ro=0 = ap pgmn = 0 (2.7) 


In the multi-centered case, we need to solve the full equations for the moduli as functions of 
the 3d coordinates x. For single-centered solutions, the moduli only depend on r; to satisfy 
the above equations, the motion of the moduli must follow null geodesics inside M3p. 


A generic null geodesic flows to the boundary of the moduli space M3p. A single- 
centered attractor flow is defined as a null geodesic that terminates at a point on the U — 
—oo boundary and in the interior region with respect to all other coordinates. This is 
guaranteed for the BPS attractor by the constraints imposed by the supersymmetry. To find 
the single-centered non-supersymmetric attractor flow, one needs to find a way to construct 
null geodesics and a constraint that pick out the ones that stop at this specific component of 
the boundary. In the next section, we will show that we can do this for models with special 
properties, and the method can be easily generalized to find the multi-centered attractor 
solution. 


2.3 Models with Symmetric Moduli Space 


The problem of finding such a constraint in a generic model is not easy. To simplify, we study 
any model whose moduli space is a symmetric homogeneous space: M3p = G/H. When 
Msp is a homogeneous space, the isometry group G acts transitively on M3p. H denotes 
the isotropy group, which is the maximal compact subgroup of G when one compactifies 
on a spatial isometry down to (1,2) space, or the analytical continuation of the maximal 
compact subgroup of G when one compactifies on the time isometry down to (0,3) space. 
The Lie algebra g has the Cartan decomposition: g = h @ k where 


(h,h])=h  {h,k] =k (2.8) 


When G is semi-simple, the homogeneous space is symmetric, and 


k,k] = h (2.9) 


The models with symmetric moduli space includes: D-dimensional gravity toroidally 
compactified to four dimensions, certain models of 4D M = 2 supergravity coupled to vector- 
multiplet, and all 4D M > 2 extended supergravity. The entropy of the last two classes is 
U-duality invariant. In the present paper, we will only consider the first two classes, namely, 
the D-dimensional gravity toroidally compactified to four dimensions, and the 4D M = 2 
supergravity coupled to ny vector-multiplet. 


The left-invariant current is 
J=M "dM = J+ Jn (2.10) 


where M is the coset representative, and J, is the projection of J onto the coset algebra k. 
The lagrangian density of the sigma-model with target space G/H is given by J, as: 


L = Tr(Jk A *3Jk) (2.11) 


The geodesic of the homogeneous space written in terms of the coset representative is 
simply 
M = M" with kek (2.12) 


where Mo parameterizes the initial point, and the z is for later convenience. A null geodesic 
has zero length: 
Ik? =0 (2.13) 


Therefore, in a homogeneous space, we can find the null geodesics that end at an attractor 
point by imposing the appropriate constraint on the null elements of the coset algebra. 


Since M is defined up to the action of the isotropy group H, in order to read off the 
moduli fields from M in an H-independent way, we construct the symmetric matrix using 
the metric signature matrix So: 


S = MSo MT (2.14) 


In all systems considered in the present paper, H is the maximal orthogonal subgroup of G 
with the correct signature: 


HSHT = S for VHEH (2.15) 


That is, the isotropy group H preserves the symmetric metric matrix Sp. Therefore, S is 
invariant under M — MH with H € H. Moreover, as the isotropy group H acts transitively 
on the space of of matrices with a given signature, the space of possible S is the same as the 
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symmetric space G/H. That is, the moduli of G/H can be combined into the symmetric 
matrix S. And the current of S is 
Js = S'dS (2.16) 


It is easy to perform the projection onto the coset algebra k. The (generalized) orthogo- 
nality condition of the isotropy group H can be expressed in terms of the subalgebra element 
h which is in H = e” as 

hSo + Soh” = 0 Vheh (2.17) 


In other words, (hSo) is anti-symmetric: (hSo)" = —(hS). Thus the coset algebra, being 
the compliment of h, can be defined as the k with (kSo) being symmetric: (kSo)7 = (ko), 
i.e. 


kT = S hsp Vkek (2.18) 
Therefore, the projection of an element g in g onto the coset algebra k is: 
Sig Sg! 
Gea 2 Oe eo 0 (2.19) 
For the left-invariant current J = M~'dM, the projection onto k is: 
J + SoJ? Sy" 
E eda (2.20) 
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It is straightforward to show that the current constructed from S' is related to the pro- 
jected left-invariant current J, by: 


Jg = S~'dS = 2(So MT)! Jk( SoM") (2.21) 
The lagrangian in terms of S is thus L = +Tr(Js A^ *3Jg). That is, the lagrangian density is 
1 
L= 128s -StV 8S) (2.22) 
which is invariant under the action of the isometry group G: 


S — GSG where GEG (2.23) 


and whose conserved current is: 


J= SVS (2.24) 


where we have dropped the subscript S in Js, since we will only be dealing with this current 
from now on. The equation of motion is the conservation of the current: 


Vel =V-(8 VS)=0 (2.25) 


We now specialize to the single-centered solutions: they correspond to geodesics in the 
coset manifold. The spherical symmetry allows the 3d metric to be parameterized as 


ds? =C(r) dr (2.26) 


Then the equations of motion involve the operator d,r?C(r)d,, and reduce to geodesic equa- 
tions in terms of a parameter 7 such that 


d 

T= rO(r). (2.27) 
dT 

The function C(r) is then determined from the equations of motion of 3d gravity. The 

equations of motion can be written as 


d aS 
“(9 )=0 2.28 
a i? ( ) 
In the extremal limit the geodesics become null, the 3d metric is flat and 


(2.29) 


In the search for multi-centered extremal solutions, where the spherical symmetry is 
absent, it is very convenient to restrict to solutions with a flat 3d metric. This is consistent 
with the equations of motion as long as the 3d energy momentum tensor is zero everywhere: 


Tx» = Tr(JaJy) = 0 (2.30) 


The coupled problem with generic non-flat 3d metric is much harder, and exact solutions 
are hard to find unless a second Killing vector is present. 


Since different values of the scalars at infinity are easily obtained by a G transformation, 
to start with, we will consider the flow starting from Mp = 1, and generalize to generic 
asymptotic moduli later. For a single-centered solution, the flow of M is M = Moe*7/?. 
Since all the coset representatives can be brought into the form e? with some g € k by an 
H-action, we can write Mo = e9/2, so M = e9/2e*7/2, And the flow of S is 


S(T) = ee" e? So (2.31) 
The charges of the solution are read from the conserved currents 


Soe 9? ke? So En 
F 
r2 


J(r) = SVS = (2.32) 


3 Torus Reduction of D-dimensional Pure Gravity 


Now we use the method introduced in the previous section to analyze pure gravity toroidally 
compactified down to four dimensions. We explain why the attractor flow generator, k, 
needs to be nilpotent, and we find the Jordan forms of k? and k. Using this information, we 
construct single-centered attractor flows. We then generalize to multicentered black holes 
in pure gravity and show that these solutions have mutually local charges and no intrinsic 
angular momentum. 


3.1 Kaluza-Klein Reduction 


The simplest example of a system that admits a 3d description in terms of a sigma model 
on a symmetric space is pure gravity in D dimensions, compactified on a D — 4 torus. The 
KK reduction to 4D parameterizes the metric as 

1 


ydet p 


Here y? are the torus coordinates, x” the coordinates on R*", and p,, the metric of the torus. 
A 4D metric with one timelike Killing spinor is then parameterized as 


dsh = Ppq(dy” + Alda") (dy? + Ajda") + ds? 1<pq<D-—4 (3.1) 


; 1 
dsj = —u(dt + widz)? + 483 (3.2) 


where u = e? , to connect with the parametrization in the later part of the paper; and 
i = 1, 2,3 denote the 3d space coordinates. 


The two expressions combine as 


ds}, = Gap(dy® + @2da*) (dy? + oda’) + ds? 0<a,b<D-4 (3.3) 


1 
— det G 
Here y? are the torus coordinates plus time, zê coordinates on R? and 
rAg 
C= ( Ppq PprA9 A ) (3.4) 


Ap rq AgPrsAg — WEET 


a 


and @* = (?, &?) is: 


OP = (AP — Abw,)dz’ D =w (3.5) 


If the forms w* are dualized to scalars a, as 


dag = — det GGap *3 da” (3.6) 


10 


the various scalars can be combined into a symmetric unimodular (D — 2) x (D — 2) matrix 


1 1 
g= ( Gab + seta det Ga ) (3.7) 
det GO det G 


In terms of the 4D fields that is 


snd es oe S 
Ppq ude p PO Ppr Ag ur/det p uv'det ap 
= r a r Se ee E ate 
S = Ao Pra uap 0a AoprsAo vdet p uydei eo udeg 0 (3.8) 
~ uydet pea ~ uydet Fel ~~ uy/det p 


The equations of motion derive from the lagrangian density L = TrVSS-!VSS"1, in- 
variant under S — U? SU for any U in SL(D — 2). As this SL(D — 2) action is transitive 
on the space of matrices with a given signature, the space of possible S is the same as 
the symmetric space SL(D — 2)/SO(D — 4,2). Notice that the signature of the stabilizer 
SO(D — 4,2) is appropriate for the reduction from (D — 1,1) to (3,0) signature. The usual 
reduction from (D — 1,1) to (2,1) would give a SL(D — 2)/SO(D — 2), while the Euclidean 
reduction from (D,0) to (3,0) gives SL(D — 2)/SO(D — 3,1) P7]. 


The coset representative under the left SO(D — 4,2) action can be described in terms of 
a set of vielbeins 


e^ = EA(dy* + w%dz') el = T ares) (3.9) 
as P 
M= ( o ) (3.10) 
det dak 
Then the symmetric SO(D — 4, 2) invariant matrix 
S = MS)M* (3.11) 


can be used to read off the solution more easily. Without loss of generality, we can take Sp 
to be the signature matrix: 


So = Diag(n, —1) = Diag(1,--- ,1,—1,—1) (3.12) 


The equations of motion are equivalent to the conservation of the SL(D — 2) currents 
J = S~'dS. Some of those currents correspond to the usual gauge currents in 4D: the first 
D — 4 elements of the last column J;,p—2 are the KK monopole currents, the first D — 4 
elements of the row before the last Jp_3, are the KK momentum currents and the element 
Jp—3,p—4 is the current for the 3d gauge field w. Regular 4D solutions must have zero sources 
for this current, otherwise w will not be single valued. 
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3.2 Nilpotency 


We now show that all the attractor flows are generated by the nilpotent generators in the 
coset algebra. To get extremal black hole solutions with a near horizon AdS x S$”, the 
function u must scale as r? as r goes to zero while the scalars go to a constant. This makes S 
diverge as 5. The most natural way for S to diverge as T° for large 7 is that k is nilpotent, 
with 

k? =0. (3.13) 


This will be the crucial condition through the whole paper. 


3.3 A Toy Example: Hyperkähler Euclidean Metrics in Four Di- 
mensions 


Not every null geodesic corresponds to extremal black hole solutions. Let’s consider a simple 
example: hyperkähler euclidean metrics in 4D. 


Although this example is not about a black hole, it is still quite instructive. The 3d sigma 
model is SL(2)/SO(1,1), i.e. AdS. The coset representative is written as 


M= ( ul? H ) (3.14) 


and the symmetric invariant 


S= ( une = ) (3.15) 


A geodesic is the exponential of a Lie algebra element in the orthogonal to the stabilizer. 
The stabilizer SO(1,1) is generated by ot. A null geodesic is hence the exponential of 
k = o? tio*. This is a nilpotent matrix, k? = 0, hence M = 1 + 7k/2. Take: 


O 3 ae 1 1 
k = o? + io = ( ie aj (3.16) 


then 


M= ( a 1 ee ) eae 


and we can read off the geodesic solution from the invariant 


ee ee) =) er T 
e a eii 
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Hence 


(3.19) 


Dualizing a, we get 

1 

+ dw = —— = —dr wt=14+- (3.20) 
r 


The 4D Euclidean metric is just the Taub-NUT metric. Notice that the multi-centered 
Taub-NUT generalization of the metric is obtained by replacing T above with some harmonic 
function >>, a The sigma model equations of motion are equivalent to the conservation 
of the current J = S~'VS, and if S is given as above with T = T(7) then the equation of 
motion are 


V?r(z) =0 (3.21) 


3.4 Single-centered Black Holes in Pure Gravity 
3.4.1 Constructing the flow generator k 


Now we look in detail at the single-centered black holes in pure gravity. Notice that as u 
goes to zero S tends asymptotically to a rank one matrix 


1 ApQg App Qp 
S = ————= | Aq aoao Qo (3.22) 


uv/det p Oy w l 


hence the matrix k? should also have rank 1. By inspection of S it is clear that a k? of rank 
higher than one gives a geodesic for which the matrix elements of p also diverge as 7° so that 
the scalar fields do not converge to fixed attractor values. 


Notice that if k is nilpotent, then S is a polynomial in Z, and the various scalars in the 
solution will all be simple functions of r for such extremal solutions! 


The explicit form for k in terms of the charges is then straightforward to write. Consider 
the Jordan form of k?: as it is nilpotent, the eigenvalues are all zero. As it is rank one, it 


has one single indecomposable block of size two: ( ; i: It is written as 
k” = —nvv" f (P,Q) (3.23) 


with v null in the metric 7, and f (P,Q) any degree-two homogeneous function of the charge 
(P,Q). This form is chosen so that v does not scale with the charge (P, Q). 
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0 1 0 
Then k must have a Jordan form with all eigenvalues zero, one block of size 3: | 0 0 1 
0 0 0 


and possibly some other extra blocks of size two. Alternatively, there is a subspace V anni- 
hilated by k, a subspace V’ whose image under k sits in V and has the same dimension as 
V’, and a single vector w such that kw € V’ and is non-zero. 


kw C V’, kV’ C V, kV =0. (3.24) 


From the symmetry of 7K, it follows that the space kV’ is made of null vectors only, and 
that kw is orthogonal to it. Because nk? = —vuT, (kw)'nkw is negative. Because of the 
signature of 7 it is straightforward to see that V’ can be of dimension at most one; hence 
there are no blocks of size 2 in the Jordan form of k. 


Taking this into account, the final form of k is simply 
k = nwu? + now" (3.25) 


where v and w are two orthogonal (D — 2)-dimensional vectors with v being null and w 
having norm —f(P, Q): 


wnv = 0, une = 0, wwe = -f (P,Q). (3.26) 


Using the fact that in k, the first D — 4 elements of the last column K; p-2 are the magnetic 
charges, and the first D — 4 elements of the row before the last Kp_3,; are the electric 
charges, and the element Kp-3,p-4 is the Taub-NUT charge, which has to vanish, we have 
2(D — 4) + 1 = 2D —7 conditions. Together with the three constraints coming from the 
norms and orthogonality condition (8.26), they can be used to solve for the 2D — 4 degrees 
of freedom in (v, w). 


The full solution of (v, w) requires one to solve some degree-four equations, hence we'll 
leave it in a slightly implicit form. Let (p,q) be two (D-—4)-dimensional vectors proportional 
to the magnetic and electric charges, so that the magnetic charge and the electric charge 
(P,Q) of the 4D gauge fields are 


P= V/p?+p-qp Q=v@+p:qq (3.27) 


And we choose f (P,Q) to be 
f(P,Q)=p-g (3.28) 


The solution of v and w written in terms of (p, q) is 


1 q+p Dq q= pP 
= syi Pe |; w= | -vč +pg |. (3.29) 
V = pP+p-q VP +p:q 
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and k can be written as 


qq — pp" VP +p-dq Vp? +p-4aP 
k= @tp-qq —-(¢+p-q) 0 (3.30) 
-yp +p- qp" 0 pet+p-q 


3.4.2 Full flow 


First, for the full flow starting from Mp = 1, the scalars for the attractor solution generated 
by this k can be read off from S(T) = e*7, by comparing with the form of S in terms of the 
AD fields: . : 

wu? = [1+ (p +p- q)(T + PA +U +p-a)(r+ 217) (3.31) 


and 
(qq® — pp™)r + [(p? +p- gad’ — lp +) + oN 


=1+ 
’ T+ (+p: a)(r + Br) 


(3.32) 
Notice that as T > 00, 7~?e~¥ has the correct limit 4“, which is the entropy where P and 
Q are the physical electric and magnetic charges. 


To generalize to arbitrary asymptotic moduli, M(r) = e9/%e*"/?, and the flow of S is 
(2.31), which can be written as S(T) = e*() Sy, where K(r) is a matrix function. From now 
on, we use lower case k to denote the coset algebra that generates the attractor flow, and 
capital K to denote the function which we exponentiate directly to produce the solution. 


We will choose K(7) to have the same properties as the generator k: 
K?(r) =0 and  K?(r) rank one (3.33) 


The equations of motion V - (S~'VS) = 0 then simplify considerably with this ansatz. If 
one further requires that the subspace image of K?(T) remains constant everywhere, such 
that 


K?(r)VK(r) = VK(r)K?(r) = 0 (3.34) 
then the current reduces to 
Jas yss (vxo) F SVK(7), KC) Sp (3.35) 
and the equations of motion are 
V?K(T) + SVK (7), K(r)| =0 (3.36) 


which is solved by a harmonic K(r). 
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It might appear hard to find a K(r) that is harmonic and satisfies all the required 
constraints. However, by remembering that the constraints dictate K(7) to have the form: 
K(r) =V WT +W VT (3.37) 


with V being null and doesn’t scale with the charge (P,Q), and W orthogonal to V ev- 
erywhere, one can simply pick a constant null vector V = v’ and a harmonic vector W (rT) 
everywhere orthogonal to v’: 


W(r)=wr+m with v-W(r)=0 (3.38) 


Here m is a (D — 2)-vector and contains the information of asymptotic moduli. Thus an 
appropriate K (r) is built: 

K(T)=kT+g (3.39) 
where 


tr IT 


ki = nu'w” +nw'v g = m'm” + nmv™ (3.40) 


Now we need to solve for (v’, w’) for the same charge (P,Q) but in the presence of m. 


The form of g guaranteed that 
[k,g] =0 (3.41) 


where we have used the fact that v’ is null and w’ is orthogonal to v’. Therefore, shifting the 
starting point of moduli does not change the current as a function of (v, w): 


k! 1 I\T Ma t\T a 
I(v',w!) = So (5 :) SR (eee) So (3.42) 
r? r 
Thus, the solution of (v’, w’) in terms of charges solved from the current does not change as 
we vary the starting point of the flow, i.e. they do not depend on the asymptotic moduli: 


v'(Q) = o(Q) w'(Q) = w(Q) (3.43) 


In summary, the flow with arbitrary starting point is simply generated by 
K(r) = nuW (rt)? +W (ryw with W =wr+m (3.44) 
where (v, w) only depend on the charges (P, Q) and m gives the asymptotic moduli. 


3.4.3 Example: 5D pure gravity compactified on a circle. 


Consider for example the case of extremal black holes in D = 5 pure gravity compactified 
on a circle. The 3d sigma model is SL(3)/SO(1,2). The symmetric invariant is 


PpP a po — T2100 =i 
Sgr = pAo — ugo p(Ao)? = Ws = Tp% =i (3.45) 
-upe ap% E 
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Then we can calculate S = Soe"? using (8.30) and compare the result to Sgr above to solve 
for all the scalars. We find that 


e™ = y/[L + (@? + pg) (r + aot + (p° + pq)(T + 5) (3.46) 


1 2 4 P42 
gya EU OE aT) (3.47) 
1+ (p? + pq)(T + 7?) 


when starting from the identity. If we allow arbitrary g the flow is too complicated to write 
explicitly here, but the attractor value of p is the same: q/p. 


3.5 Multi-centered Solutions in Pure Gravity 


In the context of pure gravity compactified on a torus, we can also give some examples of 
multi-centered solutions in the same spirit as the ones for BPS solutions in N = 2 super- 
gravity, though some important features of the latter are not present here. 


We are interested in solutions given in terms of harmonic functions which can generalize 
the single-centered extremal solutions presented above. Similar to the single-centered case, 
we exponentiate a matrix function K(Z): 


S(Z) = Xs, (3.48) 
We will choose K(Z) to have the same properties of the generator k: 


k3(z#) =0 and  K?(#) rank one (3.49) 


Using a similar argument to the single-centered flow, we require that the subspace image 
of K?(#) remains constant everywhere, such that K?(Z)VK(Z) = VK(Z)K?(#) = 0, then 
the equations of motion are 

1 
V?K(#) + zV K), K(z)| =0 (3.50) 
which is solved by a harmonic K(Z). 

A multi-centered K (7) that is harmonic and satisfies all the required constraints can then 
be built in the same way as the single-centered one: 

K(2) = nuW (2) + nW(2)v" (3.51) 


where v is the same constant null vector as in k, and W (7) is the multi-centered harmonic 
function: 


w@=> = +m (3.52) 


|Z — 
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where 7; is the position of the ith center, and w; is determined by the charges at the ith 
center, and m is related to the moduli at infinity. Requiring W (7) to be orthogonal to v 
everywhere gives the following constraints on the {w;,m}: First, taking 7 to infinity, it gives 


v:-m=0 (3.53) 


Second, w; are orthogonal to v: 
v- wi =0 (3.54) 


In addition to the constraint from the zero Taub-NUT charge condition —bz — cy = 0, this 
makes the space of each possible w; only (D — 4)-dimensional. That is, though W is a 
(D — 2)-vector, one has only (D — 4) independent harmonic functions to work with, because 
of the orthogonality to v and the requirement of no timelike NUT charges. This makes the 
solution relatively boring. 


The multi-centered solution in pure gravity does not have the characteristic features of 
the typical BPS multi-centered solution in M = 2 supergravity, where many centers with 
relatively non-local charges form bound states which carry intrinsic angular momentum. 


The basic reason is that when the ansatz K(Z) = nuW7(#) + nW(Z)v7 is used, the 
second term of the conserved currents J = So(VK + 3[VK, K])So drops out. The first 
result is that the charges of the various centers in the solution can be read off directly from 
(v,w;), and they do not depend on the positions, charges of the other centers. Thus, there 
is no constraint on the position of each center as in the MN = 2 BPS multi-centered solution; 
centers can be moved around freely. 


Moreover, the condition of no timelike Taub-NUT charge is a linear constraint on the 
charges at each center which results in a static 4D solution, as *dw = 0 leads to w = 0. 
Therefore, no angular momentum is present. 


4 Attractor Flows in Gy2)/(SL(2,R) x SL(2,R)) 


We now tackle a more complicated subject: M = 2 supergravity coupled to one vector 
multiplet. First, we reduce the theory down to three dimensions and derive the metric 
for the resulting moduli space, which is the coset Gg) /(SL(2,R) x SL(2,R)) 2 We then 
discuss the Cartan and Iwasawa decompositions of the group G'g(2), which we use to construct 
the coset algebra and translate the flow of coset representative into the flow of the moduli 
fields, respectively. We then specify the representation of G2) we will be working with, and 
describe the form of attractor flow generators in this representation. 


?Other work on this coset space has appeared recently, including [28] [29] [30]. 
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4.1 The moduli space M3p 


The 3d moduli space for N = 2, d = 4 supergravity coupled to ny vector multiplets is well- 
studied, for example in [31] [82 [33]. Some of the main results are compiled in the Appendix. 
Here we briefly review the essential points. 


The bosonic part of the action is: 
1 : = 
$= - | dav g® [r — 2ggdž ^ *4d7 — F! AGr (4.1) 
T 


where I = 0,1...ny, and Gr = (ReN )rzF? + (ImN)ry* F7. For a model endowed with a 
prepotential F'(X), 

(ImF - X);(ImF - X); 
X-ImF.X 
where Fry = 0;0;F(X). We reduce to three dimensions, dualizing the vectors (w,A/) to the 

scalars (o,B;), and renaming Aj as A’. The metric of Młp is then 


Nig = Fry + 2i (4.2) 


1 — 

d? = dU-dU+ je (do + A'dB, — BdA’) - (do + A'dB; — BAA’) + giz(z, Z)dz' - dz! 
1 

+e [UmNv -1 (dB; + Nd A*) - (AB; +N jrdA”)] (4.3) 


It is a para-quaternionic-Kahler manifold. Since the holonomy is reduced from SO(4ny + 
4)) to Sp(2,R) x Sp(2ny + 2, R), the vielbein has two indices (a, A) transforming under 
S'p(2, R) and Sp(2ny + 2,R), respectively. The para-quaternionic vielbein is the analytical 
continuation of the quaternionic vielbein computed in [34]. The explicit form is given in the 
appendix. 


For ny = 1, X! = (X°, X!). For our purpose, we choose the prepotential 
(ATP 
F(X)=- xo (4.4) 


The metric of M3, with one-modulus is with gzz = af and N and (ImN)~! being 


N= ( —(2x — iy)(x + iy) 3a(a + iy) ) m=- 1 x 


3x(x + iy) —3(2z + iy) yy \ au 3a? +y? 


The isometries of the M3, descend from the symmetries of the 4D system. The gauge 
symmetries in 4D give shifting isometries of M3,, whose associated conserved charges are: 


qrdt = Jar = Pyr — ByP;, p'dr = Jg, = Pg, + A'P;, kdr = J, = P, 
(4.5) 
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where the momenta {P,, Par, Pp,} are 


1 
P, = 5e “(do + A'dB; — B;dA’‘) (4.6) 
Pyr = e 1(ImN) 13d A? + (ReN)13(ImN~')/* (dB + (ReN) «1d A*)] — B,P{4.7) 
Pg, = eT [Im (dB; + (ReN)sxdA*)| + ATP, (4.8) 
Here 7 is the affine parameter defined as dr = — *3 sin@d0dd. (p°, pt, qı, qo) are the D6- 


D4-D2-D0 charges, and k the Taub-NUT charge. A non-zero k gives rise to closed time-like 
curves, so we will set k = 0 from now on. 


Note that the time translational invariance in 4D gives rise to the conserved current 
Jy = Py +20J, + A! Jar + Bldg: (4.9) 


where Py = 2dU. The corresponding conserved charge is the ADM mass: 2Mapydt = Jy. 


4.2 Extracting the Coordinates from the Coset Elements 


The metric for the case ny = 1 describes an eight-dimensional manifold with coordi- 
nates 6” = {u, x,y, ø, A?, A’, Bı, Bo}. This manifold is the coset space Gay/(SL(2,R) x 
SL(2,R)). The root diagram for the Cartan decomposition of G(2) is shown in Figure} The 
six roots that lie on the horizontal and vertical axes {L}, L}, LF, L3} are the six non-compact 
generators of the subgroup H = SL(2,R), x SL(2,R),: 


T Lijo = FL [Litre Lij = T 2L jv (4.10) 


> 


and the two vertical columns of eight roots {a4} are the basis of the subspace K. {a14, a24} 
for each A is a spin-1/2 doublet under the horizontal SL(2, R): 


(a(g) (aa) ae) 


And {aa1, a2; a3; Qag } for each a span a spin-3/2 representation of the vertical SL(2, R): 


Qal — 3401 Qal 0 Qal —Ga2 
[sy Aa2 = 2 Lt, Aea2 j= 341 (Zc, Aea2 |= — 2093 

Qa3 5 4a3 Ae3 2492 Qea3 —3aa4 

Gad SoA Gad Qe3 Gad 0 


All the commutators can be easily read off from the Root diagram (I), we will only write 
down the following ones which will be useful later. 


1 _ 
[a11, @14] = gla, ay3] = —4L; [a21, a24] = — 3 [az2, az] = —4L;, (4.11) 
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a21 t ari (J 
+ 
Li @ 
a22 e ai2 (J 
I 
= 3 + 
=E e Lt 
a23 © ai3 © 
Lo @ 
a24 © aig @ 


Figure 1: Root Diagram of Cartan Decomposition of G2) 


Being semisimple, the algebra of G2) has the Iwasawa decomposition g = h@a®n, where 
a is the maximal abelian subspace of k, and n is the nilpotent subspace of the positive root 
space Lt of a. In Figure P| we show the Iwasawa decomposition of Ga). The two Cartan 
generators in a are {u, y}, and {x, ø, A°, A!, B;, Bo} span a nilpotent subspace n: n” = 0 
for n € n. a and n together generate the solvable subgroup Solv of G, which act transitively 
on Msp = Ga)/SL(2,R) x SL(2,R). In particular, y generates the rescaling of y, and 
{u, x, ø, A°, A', Bi, Bo} generates the translation of {U, x, o, A°, A’, Bi, Bo}. The moduli 
space M3p can be parameterized by the solvable elements: 


Ele) = eUut(Iny)y o1x+A"A"+BrBr+00 (4.12) 


The origin of the moduli space 


a=A°=-A'=B,=By=0 z=0 y=u=l (4.13) 
correspond to X(¢) = 1. 


In Fig 2] the isometries are plotted according to their eigenvalues under the two Cartan 
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A’ o Bo © 
x © 
â o B, © 
y 
e ee ae 
6, © Ale 
Ro 
By © a’ © 
Figure 2: Root Diagram of Isometry of Msp = Gx9/(SL(2,R) x SL(2,R)). 
{u, y, x, o, A?, A!, B;, Bo} generates the solvable subgroup. 
generators u and y [BI]. {u, y} are related to aaa bl 
1 1 
u = —z[(an + a24) — (a3 + a22)] y = g[8lan + az) + (ars + a22)] (4.15) 


The three generators {o, u, } on the horizontal axis and {x, y, x} on the vertical axis form 


the horizontal and vertical SL(2, R), respectively. The vertical SL(2, 
invariance. Denote the two vertical columns of eight isometries as 


fa &n —A® Bo 
fn €2 |_| —3A' 3B, 
é €13 | B; —A! 
Ez Eis -Ê A? 


R) generate the duality 


{§14, E24 } for each A span a spin-1/2 representation of the horizontal SL(2, R), and {€a1, 02, a3, Eaa } 


for each a span a spin-3/2 representation of the vertical SL(2, R). 


3The matrix representation of u and y are 


22 


1 1 


Parameterizing the coset representative M as the solvable elements, the symmetric matrix 
S can be written in terms of the eight coordinates ġ”, from the solvable elements © 


S($) = E(P) El) (4.16) 


The coordinates are read off from the symmetric matrix S. 


4.3 Nilpotency of the Attractor Flow Generator k. 


The near-horizon geometry of the 4D attractor is AdS x S?, i.e. 


e` > \/Vey|.7 as T — œ (4.17) 
In terms of the variable u = e24 
u > ——T” as T — œ (4.18) 
BH 
The solvable element is 
M = Ut ~ u3" (4.19) 


As the flow goes to the near-horizon, u — 0 
M(t) wu? w rf (4.20) 
where —£ is the lowest eigenvalue of u. That is, M (rT) is a polynomial function of 7. 
On the other hand, since the geodesic flow is generated by k via 
M(r) = M(0)e*"”? (4.21) 


i.e., M(r) is an exponential function of r. To reconcile the two statements, k must be 
nilpotent: 
kt! =O (4.22) 


That is, the element in k that generates the attractor flow is nilpotent. Moreover, by looking 
at the weights of the fundamental representation of Ga, we see that l = 2 


k? =0 (4.23) 


Moreover, the nilpotency of the attractor flow generators guarantees that it is null: 
k’ =0 => (k? =0 => Tr(k?’) =0 (4.24) 


which means that k is null. 
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4.4 Properties of Attractor Flow 


The scalar moduli space is parameterized by a symmetric 7 x 7 matrix S which sits in Ga), 
i.e. preserves a non-degenerate three form w,j, such that nis = WijkWstuWmnot "0"? is a 
metric with signature (4,3) normalized so that 7? = 1. To facilitate the comparison with the 
pure 5D gravity case we decompose 7 as 363+ 1 of SZ(3) and pick as non-zero components 
of w the 3A3A3,3A3A3 and 3@3@1as 


1 
w = dx; A dzz A dzz + dy’ A dy? A dy? — ao A dy? A dz (4.25) 


T 


The resulting expression for 7) is 
n= dady" — dz” (4.26) 


We know that k must be an element of Ga), hence also of SO(4,3). As in the pure 
gravity case, we choose the representation such that 


SokSo = kt Sok? So = (kT (4.27) 
In this base a Go) Lie algebra element is given as 


Ai; Einj2k V" V2wi1 
k= | Irw, -A —/2v (4.28) 
— 2v1! V2w;2 0 


Here A is a traceless 3 x 3 matrix. S is a symmetric element in Gg) with signature 
{1, —1, —1,1, —1, —1, 1}, ie. S = MSMT with 

0 
1) 
0 


So = 


coos 
= 


0 
0 | = Diag(1, —1, —1, 1, —1, —1, 1) (4.29) 
1 


where m is the one for pure gravity. 


If the gauge field is turned off, then S is block diagonal 


Sy 0 0 
Slpo=| 0 Sz! 0 (4.30) 
0 0 1 


where Syr is the same as the one for pure 5D gravity. Turning on a non-zero 5D vector field 
corresponds to a more general S: 


S =e (S| poo)e™ (4.31) 
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with kı a Gy) Lie algebra matrix with w, equal to the fifth component of the gauge field, 
v? equal to the time component of the gauge field and ws equal to the scalar dual to the 
three-dimensional part of the gauge field. 


In this representation, (x,y) can be extracted from symmetric matrix S' via: 


_ S35(7) pe S33(7)Ss5(7) — S35(7)” 
S33(7) S33 (T)? 


glr) = (4.32) 


And u via: i 
tae E A (4.33) 


V S33 (T) Ss5(T) — S357)? 


The 4D gauge currents sit in J = S-1V S, where J,2(J31) is again the electric(magnetic) 
current for the KK photon, J32 the timelike NUT current, and Jya (J51) the electric(magnetic) 
current for the reduction of the 5D gauge field. 


2 
Jog = 3]; SIygp=V2IQ0 I = gia I = -V2Jp, Ja = V2JR, (4.34) 


Moreover, 

Jaz — J33 = 2Jy (4.35) 
We use Q to denote the charge matrix, where it relates to the D-brane charge {p°, pt, qi, qo} 
and the vanishing NUT charge k by 


(Qai; Qi; Qr2, Qaz) = (VZP, -V 3p", S41, Vo) Qu =-2k=0 (436) 


Since k is nilpotent: k3 = 0, 


1 
S = č" So = (1 + kr + 3K T) So (4.37) 


1 


7 T7? as T — O0. 
BH |» 


The AdS x S? near-horizon geometry of the 4D attractor dictates u = 
Therefore, the flow generator k can be obtained by 


k? = 2Vey|.(US|u—+0)So (4.38) 


Computing k? using S constructed from the solvable elements ©(¢) shows that k? is of rank 
two, its Jordan form has two blocks of size 3. It can be written as 


k? = Va; Cab S0 (4.39) 


a,b=1,2 


25 


with va null and orthogonal to each other: vg - Up = vi Sov, = 0, and Cap depends on the 
particular choice of k. Thus k can be expressed as: 


k= > (vaw? + wave) So (4.40) 
a=1,2 
where each wa is orthogonal to both vg: Wa: Up = 0, and wa satisfy 


Wa * Wh = Cab (4.41) 


Parallel to the pure gravity case, the single-centered attractor flow is constructed as 
S(T) = eX So, where we choose K (r) to have the same properties as the generator k: 


K’(T)=0 and K? (T) rank two (4.42) 
This determines K(T) = kr + g where 
k= ` [vawe + wav] So and g= ` [vam? + Mavt] So (4.43) 
a=1,2 a=1,2 


where the two 7-vectors m,’s are orthogonal to va and contain the information of asymptotic 
moduli. Using |[k, g], g] = 0, the current is reduced to 


r2 


J (4.44) 


from which we obtain va and wa in terms of the charges and the asymptotic moduli. 


5 Flow Generators in the G(2)/(SL(2, R)x SL(2, R)) Model 


We now explicitly construct the generators of single-centered attractor flows. We start with 
the BPS flow which is associated with a specific combination of the coset algebra generators 
Qaa. It can be derived from the condition of preservation of supersymmetry. We then 
construct the non-BPS attractor flow generator in analogy with the BPS one. In Section 
5.2, we write kgps and knongps in terms of the va and wa vectors. This form will be especially 
helpful in generalizing to the multi-centered case. 


5.1 Construction of Flow Generators 
5.1.1 Constructing kgps using supersymmetry 


To describe BPS trajectories it is useful to remember that the stabilizer of S in Ga) is 
SO(1,2) x SO(1, 2), corresponding to the elements of Ga) which are antisymmetric after 
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multiplication by Sọ. Geodesics are exponentials of elements that are symmetric after mul- 
tiplication by So. Such elements sit in a (2,4) representation of SO(1,2) x SO(1,2). A BPS 
trajectory is highest weight for the first SO(1, 2). Labelling the symmetric generators as aa 
under the two SO(1,2) groups, a BPS trajectory is generated by 


keps = Ce Gage A (5.1) 


The twistor z and the coefficients C4 are fixed in terms of the charges of the extremal BPS 
black hole and the condition of zero time-like NUT charge. 


To see why this is true, expand the coset element kgps that generates the BPS attractor 
flow using aaa: 
kpps = agag (5.2) 


where C®4 are conserved along the flow. On the other hand, the conserved currents in the 
homogeneous space are constructed by projecting the one-form valued Lie algebra g~! - dg 
onto k, which gives the vielbein in the symmetric space: 


g dg = GaaV™ (5.3) 
where V4 is conserved: r 
(vesd) = 0 (5.4) 


Therefore, the expansion coefficients of kgps are 
CoA = Ve (5.5) 


In terms of the vielbein, the supersymmetry condition that gives the BPS geodesics are 
written as [83]: 


V2 = 6 (5.6) 
That is: 
Veiza = 0 => CAz, = 0 (5.7) 
Define z% = ° zg, 
CSC (5.8) 


Therefore, the coset element kgpg is expanded by the coset algebra basis aaa as kpps = 
Ag aC rz”. 


Note that kgpg has five parameters (C4, z) where A = 1,...,4. As will be shown later, 

z can actually be determined by (C4) and moduli at infinity. So the geodesic generated 

by kgps is indeed a four-parameter family. It is easy to show that kgpg is null, but more 
importantly, it is nilpotent: 

kbps =0 (5.9) 


As will be shown later, kgps indeed gives the correct BPS attractor flow. 
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5.1.2 Constructing kyonpps 


To construct the non-BPS attractor flow, one needs to find an element in the coset algebra 
distinct from kpgps that satisfies: 
kNonBPS = 0 (5.10) 


The hint again comes from the BPS generator. Note that kgps = aaa P4z° can be written 
as: E E 
kppg = e "r ko pge” (5.11) 
where kù pg spans only the right four coset generators al: 
n= aa (5.12) 


That is, kgps is generated by starting with the element spanning the four generators anni- 
hilated by the horizontal SL(2) raising operator L}, then conjugating with the horizontal 
SL(2) lowering operator L,. And it is very easy to show that (k%ps)? = 0 which proves 
(kgps)* =0. 


In Gai2)/SL(2,R)’, there are two third-degree nilpotent generators in total. And since 
there are only two SL(2,R)’s inside H, a natural guess for a non-BPS solution is to look at 
vectors with fixed properties under the second SL(2, R) group. An interesting condition is to 
have positive charge under some rotation of L}, ie. an SL(2, R) rotation of X- 4—12 agi. 
Therefore, this suggests to us to start with the element spanning the four generators annihi- 
lated by the square of the vertical SL(2, R) raising operator (£7)? and then conjugate with 
the vertical SL(2, R) lowering operator Ly: 


knonBps(2) = ee KNonBPse” (5.13) 
where 


Lees = en where a,a = 1, 2. (5.14) 


And one can show that: (khongps)> = 0 which proves (kyongps)? = 0. Moreover, (knongps)? 
is rank two. 


As long as one can pick the coefficients C°4 and the twistor z that describes the SO(1, 2) 
direction to be such that the time NUT charge is zero, this generator will give nice non-BPS 
extremal black holes. All the known non-BPS solutions may be recovered this way, and 
more, as this construction gives absolute freedom to pick the charges and moduli at infinity 
for the black hole (clearly for certain values of charges and moduli the solution will crash 
into a naked singularity, but this is to be expected from comparison with the BPS case) 
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5.2 Properties of Flow Generators 
5.2.1 Properties of kgps 


We now turn to solving for va and wa in in terms of C4 and z. First, from (439) 
we know that the null space of k? is five-dimensional and the v, span the two-dimensional 
complement of this null space. For kgps = aaaC4z® the null space of (kpps)* does not 
depend on C^. Therefore, the va depend only on the twistor z = 27/21. 


Recall that we are using the basis where k has the form (4.28). From inspection of kê pg, 
we find that (v1, v2) can always be chosen to have the formý 


vı = (Vi, =m Vi, 0) v2 = (—V2, MV2, V2) (5.17) 
where 7; is a 3d metric of signature (1, —1,—1), and V1, Vz are two three-vectors with 
MeMa=0 MeBm=0 BBall (5.18) 


Since any linear combination of (v1, v2) forms a new set of (v1, v2), this means in particular 
that any v2 + cv, gives a new v2. Looking at the forms of (v1, v2), we see that Vz is defined 
up to a shifting of Vi as V2 = V2 — cV. 


An explicit computation shows that V; and V2 are given by the twistor z and u as 
(21)? + (22) iyl g zy? 
Y=] 2P- (2P Vo = —>— | Huw- zw? |. (5.19) 
2212? ztu? + z2°u! 
where the twistor u = uw is related to c by 
1 +2cz 
1 — 2cz 


The twistor representation |] of Vı and Vz are 


Vie? = 2272" Va? = zuf + Pu’ (5.24) 


p= 3 (5.20) 


4When solving for (v, w), there are some freedom on the choice of (v1, v2) and (w1, w2): firstly, a rotational 
freedom 


Ruy Rie Ris Rig 
(v1, v2) =} (v1, v2) ( Roy Roo ) and (w1, w2) = (w1, w2) ( Roi Roo ) (5.15) 
where R is orthogonal: RRT = 1 Secondly, a rescaling freedom: 
Va — TVa and Wa > = Wa (5.16) 
5With the inner product of three-vectors defined as 
Ua * Ub = vi mvp (5.21) 
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where we have used the rescaling freedom to set zu? — z?u! to be 1. Note that for the BPS 
case, the twistor u is totally arbitrary. 


Now we solve for wa. The condition that wa are orthogonal to va dictates that they have 
the form: 


Wy = (W1,7W4, 0) W2 = (W2, m W2, 0) (5.25) 

where W, and W are linearly independent, and are related to the charges by wg - Wp = Cab: 
1 1 1 

Wi $ WwW, = gau Wi d Ws = 9012 Wo d Ws = 9 22 (5.26) 


Recall that V> is defined up to a shift by Vi: Vo = VP — cVı. The consequence is that W; is 
defined up to a shift by Wz: Wı = W? +cW3. Note that the numerical factors in front of V, 
and W are opposite. Write down (W?, W2) in terms of (C4, z): 


1 (C? + 0%) + (C1+C?)z —(C?4+ 04) + (Ct + C)z 
W? = ip (C? — C4) + (C1 — C3 )z Na —(C? — C4) + (C1 — C?)z 
i 2C? + 20?z SO" + 20?z 
(5.27) 


The twistor representations of W, and W, are 


Clu2 — C2?u! C2u2?— C3u! Cle? — aA 2A 
Wı = ( C2u2 — C3u! = C8u? — Ctut ) W2= ( C272 — C321 03,2 — Otz! ) (5.28) 


Define the totally symmetric P®%”: 
Puc p12 — œ? p12 — C3 p22 — C4 (5.29) 
Then the three-vectors (W1, W2) span the four dimensional space 


(Wi, WS) aps = (Pun Pz) (5.30) 


5.2.2 Properties of kNonBgPs 


The form of va for the non-BPS case is only slightly different from the BPS case: the two 
vectors Va can be chosen to have the form: 


vi = (VizmV,0) — v2 = (Vz, —mVa, V2) (5.31) 
The twistor representation of a three-vector v = (x,y,z) is 
x+y z 
Oy = T00 + yo3 + 201 = ( 5 yai ) (5.22) 
It’s length is 
vi mv = det(o,) =a? — y? — 2” (5.23) 
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where Vi, V2 are two three-vectors satisfying the same condition as the BPS ones (5.18). 
Again, the vectors V; and V> can be written as (5.19), and the twistor representations are 
given in (5.24) with one major difference: u is no longer arbitrary, but is determined by C°4 
as: 


u2 C2 
The form of (w1, w2) are also slightly different from the BPS one (6.25) 
wi = (W1, =n: W1, 0) we = (W2,mW2, 0) (5.33) 


The (W1, W2) can be written in terms of (C'°*, z) thus: 


(CG — ed Ore on ait (Cc?) + [C107 — CC?! + (C1?)?] 
(c10? = CMG )\2 + (CG)? | + (Gee Oe ENA C1221 En (c?) 


W. : | 
1= = 
2 AUC Bis Coa _ ig One + aa] 


(5.34) 
i Zico? at (30? _ Ca _ 2C] d- [C"z j- Cc — CG 
Wz =—-=| [012 + (302 — C7)z — 207] + [012 + C2 — C74] (5.35) 
2 2102" + Qc -s Ce \2 — c”] 


2 22 4 2 
In terms of u= 4 = Sr, the twistor representation of W; and Wz are: 


WEP = uu? + (CHu? — C?ut)z%24 (5.36) 


W3? = (2%u8 +u°2f) + (C7! — CMz — 3u!) z? (5.37) 


As a consequence, the precise value of u is an extra constraint on the vectors Wa, and there 
is only a three-dimensional space of them, with (W1, W2) a linear combination of (0, V1), 
(V1, 0) and (u®uf, 2V2). 


(Wi, W2) Nongps = M(0, Vi) + n(Vi, 0) + €(u%u?, 2V2) (5.38) 


6 Single-centered Attractor Flows in G'/2)/(SL(2,R) x 
SL(2,R)) Model 


Now that we have completely characterized the generators of single-centered attractor flow, 
we can lift the geodesics to four-dimensional black hole solutions. After some calculational 
effort, we find that the BPS solution is given in terms of harmonic functions. Next, we show 
that the non-BPS case is qualitatively different, and the final solutions cannot be formulated 
so simply. 
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6.1 BPS Attractor Flow 
6.1.1 Lifting a geodesic to 4D 


We have already noted that the flow starting from generic asymptotic moduli (£o, yo) is 
generated by M(r) = e7+9)/?, with g defined before in (4.43). The matrix g has the same 
form as k. Therefore, in the BPS case, it has the expansion 


JBPS = Aa2°G" (6.1) 


where the twistor z is the same as the one in kgpg. The flow of (x,y) and u can be extracted 
from the symmetric matrix S(T) = M(r)SpM(r)? via (£32) and (433). Using k3p, = 0 
and g? = 0, S(T) = Spe*™*9 is a quadratic function of T: 


Ss5(T)Bps = p(T) + a(t) —-1 
S35(T)sps = Ya(T) + da(T) (6.2) 
S33(T)pps = €p(T)+Cp(T)-1 


where {ag(T), Ya(T), €B(7)} are quadratic functions of 7, and {8g, ôg, Çg} are linear func- 
tions of 7: 


Gat) = -2 (H H = (HVA H = 1") (i? = (By) 
Be(r) = (H?- H*)z+ (B° -— H’) 


y(r) = -$ (HTH — HPH®)(2? — 1) + 2(H?(H? + HY) — EBC + H))2) 

p(t) = - (H! + H°) — (H? + H9) 

e(t) = (HH? — (H?)?)22 + (H H* — H?H?’)z + (H°H* — (H°)}) 

Ca(r) = 2(H°z+ H’) (6.3) 


where H^ is a linear function of r defined as H4(r) = C47 + G^. 


The attractor values are reached when T — oo along the geodesic: 


P (k? So)35 ; (k? So)33(k?S0)s5 — (k?.S9)35 
= — — a ÁÁ 4 
E (k? S0)33 eee (k? So)33 oe 
and r 
T TONE PEPEE near (6.5) 


(k?.S0)33(k2.S0)55 — (k?S0) 35 


The asymptotic moduli (£o, yo) can be expressed in terms of (G4, z) by extraction from 
S= e9 So. 
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—-1.0 -0.5 0.5 1.0 1.5 2.0 


Figure 3: Sample BPS flow. The attractor point is labeled (x*,y*). The initial points of 
each flow are given by (xı = 1.5,y1 = 0.5), (£2 = 2, y2 = 4), (x3 = —0.2, y3 = 0.1), (z4 = 
-1, Y4 = 3) 
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Using this technique, one can construct all BPS single-centered black holes. The charges 
of each black hole can be read off from the current J using (4.34). One example is given 
in Figure [B| where we parametrically plot x(7) and y(r) for a BPS black hole with charges 
(p°, p', qi, Go) = (5, 2,7, —3) and attractor point (x*, y*) = (0.329787, 0.788503). g 


6.1.2 4D solution for given set of charges 


To get the solution for a specific set of charges requires more effort. In this section, we 
present the analytical result for any set of charges (p’, qr). 


The ten parameters in kgpg and g are {z u, C4, G4}, among which the twistor u is 
arbitrary, corresponding to the freedom of the shift by (vı, w2) in the definition of (v2, w1): 
(v2, W1) > (vo + cv1, W1 — CW2). The remaining true parameters are enough to parameterize 
the four D-brane charges (p’,q;) and the arbitrary asymptotic moduli (£o, yo) under the 
condition of vanishing Taub-NUT charge and fixing uọ = 1. We now solve for kgps and 
g for the given D-brane charges and (zo, yo), using the eight constraints, namely, 4 charges 
and zero Taub-NUT charge plus 3 asymptotic moduli, to fix C4 and G4, leaving the other 
twistor z unfixed!’ 


For the sake of simplicity, we will denote kgps by k for the rest of this section. Then the 
current J(Q) = 45 gives the five coupled equations: 


1 
Q = Sok + z E» gl) So (6.6) 
In order to show that the BPS flow can be expressed in terms of harmonic functions: 
A(r)=Qr+h with Q = (p, qr) and h = (h', hr) (6.7) 


we will solve g in terms of h instead of (£o, yo). The four h’s relate to the asymptotic moduli 
by 
m=L(Q>h) ypwH=yQrh) w=u(Q—h) (6.8) 


and there is one extra degree of freedom to be fixed later. 


To evaluate [k, g], we first use the commutation relation (4.11) to obtain 


A B 
lasC4, E (6.9) 
°The discriminant of the charge (5,2, 7,—3) is positive, so this is indeed a BPS solution. 


’The twistor z can be left unfixed because we will not specify the asymptotic values of the scalars with 
translational invariance, namely, ({a, A’, Br}). Fixing them can fix the twistor z. 
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where the product between C4 and G4 is defined as (C, G) = C1G4—3C?G?+30%G?—C4G!. 
Then twisting Eq with the twistor z as in (5.11) gives the commutator of k and g with 
the same twistor z: 


[k, 9] = baat Cn agpz°G?] = (C,G)O (6.10) 
where © is defined as © = —;4ye7*" Lje. On the other hand, using (£43), 
[k, g] = (vavi — vivz)So(we - my — wy - Ma) (6.11) 


© can also be written as O = (vav? — viv?) So, and we can check that (w2-m,— w),-m2) = 
1 2 


(C, G}. 


First, separate from G4 the piece which has the same dependence on (A, z) as C4 on 


(Q, 2): 
G^=GÎ+E^ with G =C^(Q—>h,z) (6.12) 


That is, g contains two pieces: 
g=g +A with = gn =Qgaz°Gf and A= aaz" E^ (6.13) 


We need to solve for E4. 


There are three constraints from (6.8). The (zo, yo) and wo are extracted from the sym- 
metric matrix S = e9So via (4.32) and (4.33). On the other hand, requiring gives 
(Xo, Yo, Uo) in terms of h: 


ee (gi So) 35 ne (9;,50)33(9,50)55 — (9450) 35 me 1 
c= - Se 2e OO as 
(g9?S0)33 (g? S0)33 (g7 So)33(97S0)55 = (97,50) 35 
(6.14) 
2 
Therefore, defining II = (e9 — ) So, I has to satisfy three constaints: 
I3 = Is = Us5 = 0 (6.15) 


in order for (6.14) to hold for arbitrary h. Using the unfixed degree of freedom in h’s to set 
(C, Gr) = 0, becomes 


1 
The zero Taub-NUT charge condition in (6.16) imposes the fourth constraint on A: the 


(3,2)-element of So(k + $[k, A])So for arbitrary k has to vanish. Combining with , we 
have 4 constraints to fix E^ to be: 


Z 
= FE? = -Et = 
1+ 2? 1+ 2? 


(6.17) 
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The remaining 4 conditions in the coupled equations (6.16) determine C4 in the BPS 
generator kgps = dgaz*C4 to be 


— qız — 3ptz? + p?z3 
E nz — 3p'2? +p 


ee 
a+) 
(1 + 2?) 
ra re ae 


The G# are then determined by Gf = agaz4C4(Q — h,z). Using the solution of C4 and 
GA, we see the product (C4,G#) is proportional to the symplectic product of (p*, qr) and 
(hi, hr): 


(C4, GA) = <Q,h> where < Q,h >= phot p'hy — qh! — qoh? (6.19) 


142 
The condition (C4, G#) = 0 is then the integrability condition on h: 
< Q, h >= peho + phy = qıh! = qoh? =0 (6.20) 


Substituting the expressions of C4 and G^ in terms of (p!, qr) into , we obtain the 
BPS attractor flow in terms of the charges (p’,q;). In particular, the pti a values are 


ot = Pot es ee Ja(P°, p*, Fs 40) (6.21) 
P RU oy ne (oon | 
where J4(p°, p', qi, qo) is the quartic Er) invariant: 
Ja(p’, p`, 11, qo) = 3(p*q1)” — 6(p"G0) (p'a1) — (pao)? — 4(p)*qo + 4p" (a1)° (6.22) 


thus J4(p°, pt, d qo) is the discriminant of charge. The attractor values match those from 
the compactification of Type II string theory on diagonal TÊ, with q) - £. The attractor 


value of u is 1 


Uw p 
= V J4(p?, pt, L, qo) 


The constraint on h from ug = 1 is then J4(h?, h!, %2, ho) = 1. 


(6.23) 


T 


Now we will show that the geodesic we constructed above indeed reproduces the attractor 
flow given by replacing charges by the corresponding harmonic functions in the attractor 
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moduli. Using the properties of A, we have proved that, in terms of k and g, the flow of 
(x,y) can be generated from the attractor value by replacing k with the harmonic function 
kT + gn: 

ar)=a(k>kr+gn) y(t) =y"(k > kr + gn) (6.24) 


Since k and g have the same twistor z, this is equivalent to replacing the C4 with the 
harmonic function C4r + G# while leaving the twistor z fixed: 


aT) = eC" — Crt + Gc z) y(r) = y*(C4 — CÎT + G, z) (6.25) 


Since C^ is linear in Q and G# linear in h, and since z drops out after plugging in the solution 
of C4 in terms of (Q, z) and G# in terms of (h, z), this proves that the flow of (£o, yo) is given 
by replacing the charges in the attractor moduli by the corresponding harmonic functions: 


tpps(T)=2pps(Q > Qr+h) — yaps(t) = Yeps(Q > QT +h) (6.26) 
The integrability condition < Q, h >= 0, in terms of H = QT +h, is 


< H,dH >= 0 (6.27) 


6.2 non-BPS Attractor Flow 
6.2.1 Lifting a geodesic to 4D 


Similar to the BPS attractor flow, the non-BPS flow is generated by M(r) = e'7+9)/?, 
and (x,y) can be extracted from the symmetric matrix S(T) by and the relevant 
elements of S(T) are given by (6.2). The only difference is that now {ap, BB, YB, OB, €B, CB} 
are changed into the non-BPS counterparts {awyp, ONB, YNB,ÔNB, ENB; Ng}, which can be 
written in terms of H°* = C7 + G and z: 


ayp(T) Z -((H”H” — H” H”) (2? = 1) oti (HPY 2? -_ 9-H’ H? F (H”)°) 
Bual(r) = (2 -—1)\(H" — H” 2) — 32H” + 22H" + H” 
yNB(T) = z2 — 1)(2z(H" H?” — H”? H”) da H?” H”) En AUS halal Z a) 
1 
dual) = 5(e(L+2)H" + 23H” — H?) — 2H” 2+ (H? — H™)) 
enB(T) = (4HY H” — H”(H”? Ta 4H”))2? _ 9H? Hx pe (H”Y 
CnB(T) = 2(H"2 + (2H? + H”™)z + H”) (6.28) 


Note that 1 = u is fixed, independent of r. The non-BPS flow written in terms of (H°*, z) 
has the same simple form as the BPS flow, i.e. the scalars are rational functions with both 
the numerator and denominator being only quadratic. This is due to the nilpotency of the 
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generator: k? = 0. Again, the attractor values are reached when T — oo, and the asymptotic 
moduli can be expressed in terms of (G°*, z) by extraction from S = e So. 


Unlike the BPS case, there are only eight parameters in kyongps and gnonpps: the two 
twistors (z,u) and (C°*, G°*) under the constraints that 


C?2 G?? 
Ce Gp 


(6.29) 


Therefore, while kgpg and ggpg can parameterize arbitrary (p’,q,;) and (Zo, yo) while leaving 
(z,u) free, all the parameters in kyongps and gnongps, including (z, u), will be fixed. 


(%3: y3) i i i i i i i l i i i i l 1 1 1 1 l 1 
—0.5 0.5 1.0 1.5 


Figure 4: Sample non-BPS flow. The attractor point is labeled (a*, y*). The initial points of 
each flow are given by: (xı = 0.539624, yı = 5.461135), (z2 = 1.67984, y2 = 0.518725), (£3 = 
—0.432811, y3 = 0.289493), (x4 = 1.28447, y4 = 1.49815), (zs = —0.499491, ys = 0.181744) 


The attractor point in terms of C'°* is 
m 7 yyg( H% = cn . _ angl Ha =. C2) en p( Ho =¥ Ce) _ Vp Hoa a Cs) 
NonBPS Ens (H%* = Ce) YNonBPS Ene ( Ho = Cea) 


(6.30) 
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with z given by 
1 
5 (—3C™ — C” + z ((2? — 3)C™ + 3z(C? + C”) + 6C”)) =0 (6.31) 


As in the BPS case, the charges of the black hole are read off from the current J using (4.34). 
We have checked that the attractor point is a non-supersymmetric critical point of the black 
hole potential Vey = |Z|? + |DZ|?: 

oVgH = 0 and DZ F 0 (6.32) 
It reproduces the results reported in the literature [9]. An example of the non-BPS attractor 
flow is shown in Figure 4] with (p°, pt, q1, qo) = (5,2,7,3) and attractor point (x*, y*) = 
(—0.323385, 0.580375). Note that J4(5, 2, 7/3,3) < 0, so this is indeed a non-BPS black hole. 
Unlike the BPS attractor flow, all the non-BPS flows starting from different asymptotic 
moduli have the same tangent direction at the attractor point. The mass matrix of the black- 
hole potential at a BPS critical point has two identical eigenvalues, whereas the eigenvalues 
at a non-BPS critical point are different. The common tangent direction for the non-BPS 
flows corresponds to the eigenvector associated with the smaller mass. 


6.2.2 4D solution for given set of charges 


We now discuss how to construct the non-BPS black hole solution for a specific set of charges 
(p', qr) i 


One major difference between the non-BPS case and the BPS case is that 


[E NonBPS, JNonBPS] =0 (6.33) 

automatically, since the forms of (w1, w2) and (m1, M2) guarantee that w1 -M = wz: Mı = 0. 
Thus the charge equation becomes simply 

QnonBpes = SolkNonBPs) So (6.34) 


These five coupled equations determine the two twistors (z, u) and C in terms of (p7, qr). 
Similar to the BPS case, the four equations which determine the D-brane charge allow us to 
write C°% in terms of the charges (p°, p', q1, qo) and the twistor z via 


(—2qo + 6(p" — qo)z? + 4(p® + a1) z? — 6p'z*) 


œH 
VEE 

g2 (+) = 2(2pt ~ A E + Ophea? 
J2(1 + 22)? 

Ga a Wate Ae ant a 


V2(1+27)8 
o2 et (p +59) — 2(2p" — q)? = (P + BH (6.35) 
V2(1 + 27)? 
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and u = ce In contrast to the BPS case, the G°* do not enter the equations and therefore 
cannot be used to eliminate the twistor z. Requiring the Taub-NUT charge to vanish gives 


the following degree-six equation for the z: 
pz? + 6p'2° — (3p? + 4q,)2* —4(Sp" — 2q0) + (3p? + 4q1)2? + 6p'z — p° = 0 (6.36) 


The three parameters in gyongps, namely, G°* with the constraint a = u are then fixed 
by the given asymptotic moduli (xo, yo) and uo = 1. 


Similar to the BPS flow, the full non-BPS flow can be generated from the attractor 
value by replacing C°* with the harmonic function H°*(r) = Cr + G°, while leaving 
z unchanged as in (6.25). However, there are two major differences. First, the harmonic 
functions H°* have to satisfy the constraint 

H??(r) c2 G2 


-< SS YH S 


Note that this does not impose any constraint on the allowed asymptotic moduli since there 
are still three degrees of freedom in G°* to account for (£o, yo, Uo). We will see later that it 
instead imposes a stringent constraint on the allowed D-brane charges in the multi-centered 
non-BPS solution. 


Secondly, unlike the BPS flow, replacing C°® in the attractor moduli by the harmonic 
function H°*(r) is not equivalent to replacing the charges Q with H = Qr + h as in (6.26). 
The twistor z here is no longer free, but is determined in terms of the charges as a root of the 
degree-six equation (6.36), so replacing Q by QT + h, for generic Q and h, would not leave 
z invariant. Therefore, the generic non-BPS flow cannot be given by the naive harmonic 
function procedure, as proposed by Kallosh et al [19]. Next, we will define the subset of the 
NonBPS single-centered flow that can be constructed by the harmonic function procedure. 


When the attractor has only D4 — DO charges, namely, Qao = (0, pt, 0, qo), has a 
root z = 0, which is independent of the value of charges. If the asymptotic moduli h is also 
of the form of h4o = (0, ht, 0, ho), replacing Q40 by QaoT + hag would leave the solution z = 0 
invariant. Now we will use the duality symmetry to extend the subset to a generic charge 
system with restricted asymptotic moduli. 


The one-modulus system can be considered as the STU attractor with the three moduli 
(S,T,U) identified. Since the STU model has SL(2,Z)? duality symmetry at the level of 
the equations of motion, the one-modulus system has an SL(2,Z) duality symmetry coming 
from identifying the three SL(2, Z) symmetries of the STU model. That is, the one-modulus 
system is invariant under the following element of SL(2, Z)’ 


a a b a b a b ; 
Sa a 2)9(% a with ad —be=1 (6.38) 
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The modulus z = x + iy transforms as 


7 az +b 
Tz = 6.39 
a cz+d ( ) 


and the transformation on the charges is given by [35]. A generic charge (p°, p1, q1, qo) can 
be reached by applying the transformation I on a D4 — DO system. 


Under the aforementioned transformation, a D4 — D0 system transforms into ÎQ40 


0 —c(3d?p' + qo) 
1 1 2 
|p a E d(2bc + ad)p* + ac*qo 
Qu = C S PQs = 3(b(bc + 2ad)p! + a?cqo) 
qı a(3b*p! + a*qo) 


The solution of the twistor z with the new charges igi is 


+ 2 2 
pa eee (6.40) 
C 


independent of the D4 — DO charges we started with. Now given an arbitrary charge Q, 
there exists a transformation Îo such that Q = ÎoQao for some Qo. The twistor z remains 
invariant under Q — Qr + Taha for arbitrary h4o. We conclude that the non-BPS single- 
centered black holes that can be constructed via the naive harmonic function procedure are 
only those with (Q, h) being the image of a single transformation Î on the (Quo, h40) from a 
D4 — DO system: 


LnB(T) = EAT, =; TQaor + Thao) YnB(T) = Ynp(EQao = ÎQuoT T hao) (6.41) 


Since we are considering arbitrary charge system, the constraint is on the allowed values of 
h. 


7 Multi-centered Attractor Flows in G(/(SL(2, R) x 
SL(2,R)) model 


As proven in the pure gravity system, the multi-centered attractor solutions are given by 
exponentiating the matrix harmonic function K (7): 


S(#) = eX, (7.1) 
with A(Z) having the same properties as the generator k: 


(a) =0 and K?(Z) rank two (7.2) 
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We now describe how to formulate A(#) for multi-centered solutions in G'(2 


The K (2) satisfying all the above constraints is constructed as 
K (2) = X [vaWa(2)" + Walz] So (7.3) 
a=1,2 


with v, being the same two constant null vectors in k, and the multi-centered harmonic 
function 


(Wali 
W,(Z) = > Ma 7.4 
@=Dea (7.4) 
is everywhere orthogonal to va. The two 7-vectors (m1, M2) contain the torma iot of 
asymptotic moduli and has the same form as (w1, w2). Write K(z) as K(Z) = $; aA mtg 
where 


k= ` [Wal Wa) + (Wwa)iv?]So and g= ` [vam] + mav? ] So (7.5) 
a=1,2 a=1,2 
Since v only depends on the twistor (z, u), and (w1, w2) are linear in C4 or C°*, the above 
generating procedure is equivalent to replace C4 or C°% by the multi-centered harmonic 
functions while keeping the twistor (z, u) fixed. 


Next we discuss the properties of the BPS multi-centered attractor solution and non-BPS 
ones separately, since they are very different in character. 


7.1 BPS Multi-centered Solutions 


In constrast with ne multi-centered solutions in pure gravity, now the second term of the 
current J = VK + 3[VK, K] does not vanish automatically since 


eA an ea (7.6) 


Therefore, the centers are no longer free, and we cannot simply read off the charges from J. 
Instead, we need to solve for Cf and G4 in a set of 5N coupled equations. The divergence 
of the current is 


lki, kj 
V. J= a (Z — £;)So(ki + = “lk a] +5 ena 


So tt 
z zp (7.7) 
Using Q; to denote the charge matrix which relates to the D-brane charge {p°, p1, q1, qo}i as 
in (4.36), and with Q32 as the vanishing NUT charge, we have 5N coupled equations from 
Q; = = J, Vo: 
lki, kj] 


Hs 


Q: = Salhi + 5 kog += SLA (7.8) 
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The generators of the multi-centered BPS attractor solution {k;} and g have 4(N +1)+2 
parameters in total: the two twistors (z,u) and {CA,G4}. On the left hand side of (7.8), 
there are also 3N — 3 degrees of freedom from the position of the centers ṣi. On the other 
hand, a generic N-centered attractor solution has 4N D-brane charges (p’,q;), and three 
additional constraints from the asymptotic moduli (£o, yo) and uo = 1. As we will show, like 
the single-centered BPS solution, the three asymptotic moduli, together with the vanishing 
of the total Taub-NUT charge, determine the 4 G4 inside g. Moreover, as in the single- 
centered case, we can solve CA in terms of the 4D D-brane charges Q; while leaving (z, u) 
unfixed. The remaining N — 1 zero Taub-NUT charge conditions at each center will impose 
N — 1 constraints on the distances between the N-centers Ti. 


First, integrating over the circle at the infinity, }>,Q; = = fV -J gives the sum of the 
above N matrix equations: 


Qia = Y Qi = SoC + FID ka g) So (7.9) 


which is the same as the one for the single-center attractor with charge Q;.. This determines 
g to be g = gn + A, same as the one for single-centered attractor in (6.13). As in the single 
centered case, h is fixed by the asymptotic moduli (xo, yo) by 


To = Tpps(Q — h) Yo = Ypps(Q — h) (7.10) 
and the two constraints: 
< Qiot, h >= 0 Fiche h', E ho) =1 (7.11) 


We have used the vanishing of the total Taub-NUT charge to determine A. Next, we will 
use the remaining coupled 5N — 1 equations to solve for the 4N {CA} and impose N — 1 
constraints on the relative positions between the N centers {Z;} where i = 1,--- , N. 


The tentative solutions of CA are given by (6.18) with (p’, qr) replaced by (p}, qr). The 
flow generator of each center k; is then ki = aa pa Substituting the solution of k; and 
g = gn + A into (28), and using 


where all the k;’s and gp have the same value for the twistor z, we get 


< Qi, Qj > 


|Z; — £5 


Qi = Solki+ < Qh >O +- alk A] Es 3 ©) So (7.13) 


Just as in the single-centered case, the solution of k; and the form of A guarantee that 
Q; = Solki + = = [kis A) So (7.14) 
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We see that as long as the following integrability condition is satisfied: 


< Qi Qj > 


|Z; — Tj] 


<Qih>+>- =0 (7.15) 
j 


the k; and g given above indeed produce the correct multi-centered attractor solution. Just 
like in the single-centered case, the multi-centered solution flows to the correct attractor 
moduli (xž,yž) near each center, independent of the value of z. It also follows that the 
multi-centered solution can be generated by replacing the charges inside the attractor value 
by the multi-centered harmonic function: 


zsrs(ë) = shps(Q >) r gt yars(#) = ybeslQ > Daa th) (10) 


z-z |x — z;| 


The sum of the N equations in the integrability condition (2.15) reproduces the constraint 
on h: < Qio, h >= 0. Thus the remaining N — 1 equations impose N — 1 constraints on the 
relative positions between the N centers {7;} with i = 1,---,N. From and (4.6), we 
see that «dw is given by J23. Defining the angular momentum J by 


k 
Wi = 2eijn I? — as r — 00 (7.17) 
r 


we see that there exists a nonzero angular momentum given by 


Sy 


= >» as HQ Q;) (7.18) 


i<j 


Thus we have shown that our multi-centered BPS attractor solution reproduces the one 


found in [18]. 


7.2 non-BPS Multi-centered Solutions. 


For given (z,u) and raa. the non-BPS multi-centered solution is the same as the 
single-centered i asin with H°*(r) replaced by the multi-centered harmonic function 
HG) = oe paj t o" ama the constraint 


H? ( #) Ce G?? 


“HE@ CP” GB 


(7.19) 


Accordingly, the attractor values at each center is the same as (6.30) with the corresponding 
Ce*, The asymptotic moduli are obtained by extraction from = = 2 So. 
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The equation of motion for the non-BPS multi-centered solution simplifies a great deal 
since 


[Vk (z), K(z)] =0 (7.20) 
automatically, following from the fact that for the non-BPS system: 
(w1); : M2 = (w2)i: Mı = 0 (7.21) 


which are guaranteed by the forms of NonBPS (w1, wə); and (m1, M2). Therefore, the 5N 
equations decouple into N sets of 5 coupled equations: 


Qi = So(ki) So (7.22) 


Equation differs greatly from the BPS counterpart (7.8). Firstly, the generators of 
the multi-centered non-BPS attractor solution {k;} and g have 3(N +1) +2 parameters: the 
two twistors (z,u) and {C?*, G°*} with the constraint (7.19). In constrast to the BPS case, 
g does not enter the equation. Thus we can simply use the three asymptotic moduli, without 
invoking the zero Taub-NUT condition, to determine the 3 G°* inside g. Secondly, unlike 
the BPS multi-centered solution, the position of the centers 7; do not appear in the equation, 
therefore there will be no constraint imposed on them: the centers are free. Last but not least, 
the remaining 3N +2 parameters in (z, u) and C’“ are not enough to parameterize a generic 
N-centered attractor solution, which has 4N D-brane charges (p/,q7;;). Accordingly, the 
multi-centered non-BPS attractor generated by this ansatz will not have arbitrary charges. 
Combining with the fact that 7; do not appear in the R.H.S of the equation, we find that 
all the N vanishing Taub-NUT charge conditions can only act on the charges on the L.H.S. 
We conclude that, in total, there will be 2N — 2 constraints on the allowed charges. 


Now we will show in detail the derivation of the constraints. First, like in the single- 
centered NonBPS solution, the absence of the Taub-NUT charge at infinity fixes z via 


2 Qi = Soe ki) So (7.23) 


The solution is the same as the solution to with the charges replaced by the total 
charges of N centers: z = z(Q — )°,Q;). Since all the N centers share the same twistor 
z, the absence of the NUT charge at each center imposes N — 1 constraints on the allowed 
charges Q;: all z(Q;) have to be equal. 


The remaining 4N equations in determine C in terms of z and Q;. Since the 
N-centers decouple, (7.22) for each center is the same as the single-center one (6.34). Thus 
the solution of C° is given by (6.35) with (p!, qr) replaced by (p/,q7,:). Again, since all the 
centers share the same twistor u, the condition imposes another N — 1 constraints on 
the allowed charges. Solving these 2N — 2 constraints, we see all the charges {Q;} are the 
image of a single transformation Î on a multi-centered D4 — DO system Qao: 


Qi = Qu (7.24) 
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The charges at different centers are all mutually local 


(Qi, Qj) = 0 (7.25) 


Except for the constraint on the charges, the N centers are independent, and there is no 
constraint on the position of the centers. A related fact is that the the angular momentum 
is zero. 


Like the non-BPS single-centered case, though the multi-centered solution can be gener- 
ated from the attractor value by replacing e N the multi-centered harmonic e 
H(z) = DE z z] C% under the constraint ( , while leaving z unchanged as in 
the generic T cannot be generated via m S function procedure used in me 
BPS case, namely, by replacing the charges inside the attractor value by the corresponding 
multi-centered harmonic functions. The reason is again due to the fact that the twistor z, 
being a function of charges, does not remain invariant under this substitution of charges by 
harmonic functions. The multi-centered non-BPS solutions that can be generated by the 
harmonic function procedure are those with {Q;,h} being the image of a single Î on the 
{Qaoi,; hao} of a pure D4 — DO system: 


ÎQa, i 
|Z — Z| 


typ(£) = t% p(Q > oe +Thao) = -ywa(Z) = yup Qao > + Îhzo) 
Iz 


(7.26) 
It appears that the existence of a simple linear ansatz for “superimposing” single center 
solutions exists in general only for mutually local extremal black holes, and only in the 
supersymmetric case does it extend to mutually non-local centers. 


To summarize, the non-BPS multi-centered solution is different from the BPS case be- 
cause it imposes no constraints on the position of the centers, but instead on the allowed 
charges Q;: the choice of charges at each center are restricted to a three-dimensional sub- 
space, and they are mutually local. The result is that the centers can move freely, and there 
is no angular momentum in the system. It does not have interesting moduli spaces of centers 
with mutually non-local charges, so it is as “boring” as the pure gravity case. 


8 Conclusion and Discussion 


In this paper, we find exact single-centered and multi-centered black hole solutions in theories 
of gravity which have a symmetric 3D moduli space. The BPS and extremal non-BPS single- 
centered solutions correspond to certain geodesics in the moduli space. We construct these 
geodesics by exponentiating different types of nilpotent elements in the coset algebra. Using 
the Jordan form of these nilpotent elements, we are able to write them down in closed explicit 
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form. Furthermore, we can use a symmetric matrix parametrization to recover the metric 
and full flow of the scalars in four dimensions. 


We have also generalized the geodesics to find solutions for non-BPS and BPS multi- 
centered black holes. The BPS multi-centered solution reproduces the known solution of 
Bates and Denef. Given our assumption that the 3D spatial slice is flat, we find that a non- 
BPS multi-centered black hole is very different from its BPS counterpart. It is constrained 
to have mutually local charges at all of its centers and therefore carries no intrinsic angular 
momentum. It is possible that if we dropped this assumption, we could find more general 
non-BPS multi-centered solutions. Such configurations would probably be amenable to exact 
analysis only in the axially symmetric case, using inverse scattering methods. 


There are many other avenues for future work. One could explore nilpotent elements in 
other symmetric spaces, and see whether non-BPS bound states with nonlocal charges exist. 
In particular, it would be interesting to study Egg)/SO*(16), which is the 3d moduli space 
for d = 4,N = 8 supergravity. We would also like to find a way to modify our method 
so that we can apply it to non-symmetric homogeneous spaces, and eventually to generic 
moduli spaces. We could then study the much larger class of non-BPS extremal black holes 
in generic NV = 2 supergravities. 
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A Derivation of the moduli space M3p 


Here we briefly review the derivation of the 3d moduli space M3p from the c*-map of the 
4D supergravity coupled to ny vector multiplets [33]. 


The bosonic part of the action for the M = 2 supergravity coupled to ny vector-multiplets 


1 . = 
S=- | d'zyg® [R — 2ggdž ^ »,dz — F! A Gi (A.1) 
T 
where the ranges of the indices are i,j = 1,...,ny and J = 0,1,...,ny, and Gr = 
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(ReN),7F7 + (ImN) 7 * F7. The complex symmetric matrix V7, is defined by 
Fy = Niz X” D;Fı = N1yD;X? (A.2) 
For model endowed with a prepotential F(X), 


(ImF - X);(ImF - X); 


Nis = Fis +2 e F.X 


(A.3) 
where Fry = rð F(X). 


After reduction on the time-like isometry, the action is S = —¢ f dt f dx £. The 3D 
lagrangian £ has three parts: L = Loravity + Lmoduli + Lem where 
1 l aw 
Lorata = -3v8 R + dU ^ «dU — 1° dw A *dw 
Lmoduli = ggd A *dz! (A.4) 
1 1 
Lem. = 50 UmN) 150.45 A «d Ad + ze (Um) 15(d A + Ajdw) A *(dA’ + Af dw) 
+(ReN) 7A} A^ (dA? + Af dw) 
The dual scalars for w and A’ are defined by: 


eU (ImN) 17 * (dA? + Agdw) + (ReN);jdA3 = —ddar 
eV x dw + (Aida: —¢,41d45) = —dd,, (A.5) 


After renaming the variables ¢,, — c, Aj — Al, ġar — By, we obtain the 3d lagrangian in 
terms of scalars only: 


L 


1 1 
-3v8 R+dU A dU + je (do + A'dB, — BdA!) A *(do + A'dB, — BrdA!) 


, - 1 
+g(z, Z)dz’ A *d7 + 56 (Lm) 150A! ^ *d AZ 


+56 (Im) (AB; + (ReN) 1d A*) A *(dB; + (ReN) 71d A") 


1 
= y8 R+ Imn0ah" 0p” (A.6) 
where, as before, 6” are the 4(ny + 1) moduli fields: ¢” = {U, z', 2,0, Al, By}, and g,, is 
the space time metric, gmn is the moduli space metric. Therefore, the moduli space M3p 
has metric: 
1 : = 
ds? = dU -dU + A (do + A’dB; — B;dA’) - (do + A'dB; — Byd A’) + 9,;(z, Z)dz' - d7’ 
1 
+568 [ImN (dB; + Nd A*) - (dB; +N y,dA*)| (A.7) 
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It is a para-quaternionic-Kaéhler manifold. Since the holonomy is reduced from SO(4ny + 
4)) to Sp(2,R) x Sp(2ny + 2, R), the vielbein has two indices (a, A) transforming under 
S'p(2, R) and Sp(2ny + 2,R), respectively. The para-quaternionic vielbein is the analytical 
continuation of the quaternionic vielbein computed in [34]: 


iu Vv 

yes _ EP f iE* 
ih æ 
-ð iü 


The 1-forms are defined as 
e£/-U X! (a By + N71 dA") 
r efdz’ 


e VetgeX/? Ds X! (AB; + N70 A’) 


E 
MEO O I 


» = -dU+ =e (da + AldB, — BdA’) (A.8) 


where e? is the veilbein of the 4D moduli space, and the bar denotes complex conjugate. 
The line element is related to the vielbein by 


ds? = —u- a+ gge" -P — gg EB? B® + 0-0 = captan V & VB (A.9) 


where cag and cag are the anti-symmetric tensors invariant under Sp(2,R) = SL(2, R) and 
Sp(2n, + 2, R). 


The isometries of the M3, descends from the symmetry of the 4D system. The gauge 
symmetries in 4D gives the shifting isometries of M3): 


AP — A +AA! 
Br — Br+ABr (A.10) 
g- == o+ Ao PABA — AA! B; 


The conserved currents and charges are given by (4.5) and the discussion thereafter. 
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